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ABSTRACT 

Spectral observations of the massive colliding wind binary Eta Carinae show phase- 
dependent variations, in intensity and velocity, of numerous helium emission and absorption 
lines throughout the entire 5.54-year orbit. Approaching periastron, the 3D structure of the 
wind-wind interaction region (WWIR) gets highly distorted due to the eccentric (e ~ 0.9) bi¬ 
nary orbit. The secondary star (t]b) at these phases is located deep within the primary’s dense 
wind photosphere. The combination of these effects is thought to be the cause of the particu¬ 
larly interesting features observed in the helium lines at periastron. We perform 3D radiative 
transfer simulations of rj Car’s interacting winds at periastron. Using the simplex radiative 
transfer algorithm, we post-process output from 3D smoothed particle hydrodynamic simu¬ 
lations of the inner 150 au of the rj Car system for two different primary star mass-loss rates 
(Mjjjs,)- Using previous results from simulations at apastron as a guide for the initial conditions, 
we compute 3D helium ionization maps. We find that, for higher Mtja^ He^^-ionizing pho¬ 
tons are not able to penetrate into the pre-shock primary wind. He^ due to tjb is only present in 
a thin layer along the leading arm of the WWIR and in a small region close to the stars. Low¬ 
ering allows 77 b’ s ionizing photons to reach the expanding unshocked secondary wind 
on the apastron side of the system, and create a low fraction of He^ in the pre-shock primary 
wind. With apastron on our side of the system, our results are qualitatively consistent with the 
observed variations in strength and radial velocity of r] Car’s helium emission and absorption 
lines, which helps better constrain the regions where these lines arise. 

Key words: radiative transfer - hydrodynamics - binaries: close - stars: individual: Eta 
Carinae - stars: mass-loss - stars: winds, outflows 


1 INTRODUCTION 


The massive colliding wind binary Eta Carinae (t; Car) is 
one of the most luminous stellar objects in the Galaxy 
(f-Total ^ 5 X 10^ Lq). During its highly eccentric {e ~ 0.9), 
5.54 yr orbit, the slow but dense wind of th e primary star 
(Voo ai 420 kms~' « 8.5 x IQ-^ Moyr"' 


n e primary siar t/a 
: iHillier et al.ll200ll : 


Groh et al. 2012ah collide s with the faster (v op ^ 3000 kms“^; 
Pittard & Corcoranl l2002l : IParkin et ^ 120091) . but less dense 
~ 1.4x 10“^ Moyr“^) wind of the secondary star t/b- 


Every periastron, orbital motion and the binary’s large eccen¬ 
tricity highly distort the stellar winds and wind-wind interac¬ 
tion region (WWIR), causin g t/r to become deeply embedded 
within the denser t/a wind ( Okazaki et al.l |200^ : iMadural I 2 OIOI : 


E-mail: clementel@strw.leidenuniv.nl 


IParkin et ^l201ll : iMadura et~^l2012bl . l2013h . At periastron, the 
dense distorted t/a wind and WWIR trap the ionizing radiation 
from t/b that is responsible for the formation of numerous high- 
ionization (ionization potential > 13.6 eV) emission and absorp¬ 
tion line s observed during the broad pyt of the orbit around 
apastron (yerner. B ruhweile r & Gulll12 005 i lPan iineli et al.l 2008al: 

]l2009l:lMa ' 


Gull et al.ll2009t IMadurall20 id: I 

2OI2I : Madura et al.ll2012hLl2013 : lTeodoro et al.ll201.3h . 


Gull et al II2OIiL lMadura & Grolj 


In a recent paper, IClementel et al.l (l2014bL hereafter C14b) 
presented results from 3D radiative transfer (RT) simulations focus¬ 
ing on the ionization structure of helium, due to t/b ’s ionizing radia¬ 
tion, in the inner ~ 155 au o f the rj Car system at an orbital phase of 
apastron. As summarized in lC14bl. helium spectral features provide 
important information on both the geometry and physical proper¬ 
ties of the 77 Car binary and the individual stars. While various he¬ 
lium features are present throughout the entire orbit, they show their 
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most interesting behavior during periastron passa ge (iNielsen et alJ 
l2007l : lDamineli et aLll2008al : [Teodoro et alJl2012l) . 

The goal of this paper is to extend the work of lC14bl and ana¬ 
lyze the effects of t/b’s ioniz ing rad iation during the spectroscopic 
low state at periastron. As in lC14bL we focus on the inner ~ 155 au 
of the 77 Car system and apply radiative transfer post-processing 
simulations to 3D smoothed particl e hydrodynamics (S PH) sim¬ 
ulations of 77 Car’s colliding winds (iMadura et al .112013l hereafter 
Ml3) to compute 3D maps of the ionization structure of He. This 
work sets the stage for future efforts to compute synthetic He line 
profiles for direct comparison with available observations. The goal 
is to constrain where and how the observed broad He emission and 
absorption lines form during the periastron event. We further aim 
to explain the temporal behavior of the broad He lines as the 77 Car 
system moves from apastron through periastron. Again, we note 
that we focus solely on interpreting the broad emission and absorp¬ 
tion features of He that arise in the stellar winds and WWIRs, and 
not the narrower (< 50 kms“^) features that form in the Weigelt 
blobs and other dense, slow-moving near-equatorial circum stellar 
ejecta (IWeigelt & Ebersbergejl 19861 : IPamineli et al.ll 200 ^ . 

This paper is organized as follows. In Section |2l we summa¬ 
rize our numerical approach. Section [3] describes the results. A dis¬ 
cussion of the results and their implications is in Section a Sec¬ 
tion |5] summarizes our conclusions and outlines the direction of 
future work. 


2 METHODS 


The numerical approaches in this pa per, bo th for the SPH and RT 
simulations, are identical to those inlC14bl. As in prev i ous works 


us e the simplex algorithm ( Rjtzerveld^Jcke 200^: iRitzervelc 


2007 : IPaardekooper. I^uip & Ickdboiol : iKruip et ^l 2010 l : iKruip 


2011 ) to post-process 3D SPH si mulation output jM13h. F or fur- 
ther details, we refer the reader to lClementel et al.l j2014ah . lC14bL 
an d references therein. We also r etain the na ming convention used 
in lM13l.lciementel et al.l (l2014ah . and lC14bl when referring to the 
SPH and simplex simulations in this paper, namely. Case A (Mtja = 
8.5x10-4 Moyr-i), Case B = 4.8x10-4 Moyr-^), and 
Case C (MrjA - 2.4x 10-4 Moyr-^). In the following, we briefly 
describe only the key aspects and differences that are directly rele¬ 
vant to this work. 

Due to the high eccentricity of the 77 Car system (e ~ 0.9), as 
the stars approach periastron, both the 3D structure and dynamical 
evolution of the individual stellar winds and WWIR change dras¬ 
tically, compared to the majority of the orbital period. The stellar 
separation drops to ~ 1.5 au and orbital speeds increase greatly, 
becoming comparable to the wind speed of rj^ (see Section [3] for 
details). The rapid orbital motion during periastron thus leads to 
major changes in the structure of the winds and WWIR on rela¬ 
tively short time-scales of the order of a few days jOkazaki et al.l 
120081 IParkin et ^|201 ll. iMadura et al.ll 2012 bl : lM13h . 


To properly resolve this complex si tuatio n, we use two differ¬ 
ent SPH simulation domain siz es fromlM13l . The larger domain 
is the same as the one used in IC14bl (r = 10a = 155 au, where 
a = 15.45 au is the length of the orbital semimajor axis), and al¬ 
lows for a dir ect com parison of the results in this paper with those 
at apastron in lC14bl . The smaller domain simulations have a com¬ 
putational domain radius r = 1.5 a = 23 au and eight times the num¬ 
ber of SPH particles used in the r = 10a simulations, leading to a 
roughly factor of two improvement in the overall spatial resolution 
compared to the r = 10a simulations. The high-resolution simula¬ 


tions focus only on phases around periastron (0 = 0.97-1.03, where 
by convention, 0 = t/2024 = 0,1,2,..., is defined as periastron, and 
0 = 0.5,1.5,..., is defined as apastron), with output written every 
A0 8 X 10-4, which is approximately every 1.6 days of the 2024- 
day orbit. These smaller, higher-resolution simulations ensure that 
we are adequately resolving the stellar winds and WWIR directly 
between the stars, as well as any possible W WIR ‘c ollapse’ that 
may affect the escape of 773 ’s ionizing radiation (lM13h . The higher 
resolution r = 1.5 a simulations are also useful for checking the ac¬ 
curacy of the r = 10 a results. 

We consider the ionization of hydrogen and helium atoms 
by both collisional- and photoionization, a nd assume the sam e 
abundance by number of He relative to H as iHillier et HI (l200lh . 
^He/^H = 0-2. We emplo y a sin gle photoionizing source located at 
the position of 77 B. As in IC14bL we consider 77 B to be an 05 giant 
with Tgff ~ 40,000 K, we a ssume a total ionizing flux for H an d He 
of 3.58 X 104^ photons s-^ (iMartins. Schaerer & Hilliejl2005h . and 
we use three bins to sample the stellar spectrum, which we approxi¬ 
mate with a black body. The bin edges of each frequency bin are set 
by the ionization energy of each species. The first bin ranges from 
the ionization frequency of H^'*' (yH 0 + = 3. 28 X 10^^ Hz) to that 
of He^^ “ 5.93 X 10^^ Hz), the second from yjjg 0 + to yHe+ 

(1.31 X 10^^ Hz), and the third from yHg+ to a maximum frequency 
equal to ten times yH 0 +. We use an effective cross-section repre¬ 
sentation to determine the correct number of absorptions within 
each frequency bin. We set the black body temperature to the 
value that produces the correct ratio photonsH 0 +/photonsp^g 0 + (in 
th is case T^h = 49,000 K), in accorda nce with the qiio+1q^^o+ ratio 
in iMartins. Schaerer & Hilli^ (l2005h . effectively zero for 

In figures displaying our simplex results for density and tem¬ 
perature, we visualize the average of the four vertices that com¬ 
pose each Delaunay cell. Unfortunately, this approach leads to cell 
values that are difficult to interpret if the vertex values differ by 
several orders of magnitude. This might occur for the fractions of 
He^^, He^, and He^^. In order to visualize results for these quan¬ 
tities in an understandable way that more truthfully represents the 
physics of our RT simulatio ns, we adopt the visualization approach 
described in section 2.2.5 of lC14bl . In all figures showing ionization 
structure, we display the minimum vertex value for the fraction of 
He^^ and the maximum vertex value for the fraction of He^^. For 
He^, we show the maximum value whenever the temperature of the 
gas is lower than 10^ K, and the minimum otherwise (for temper¬ 
atures > 10^ K, the fraction of He^, due to collisional ionization, 
is ^ I Q-^ or lower). For the reasons discussed in section 2.3.1 of 
IC14bl (i.e. H is already ionized in the inner ~ 120-150 au region 
centred around 77 a), we do not include plots of the ionization struc¬ 
ture of H. We emphasize that the visualizations are merely to help 
guide the reader, and neith er the p hysics nor the conclusions of our 
work depend on them fsee lC14bl for details). 

Using the method described in section 2.3.1 of lC14bL we ac¬ 
count for the inner He^ volume around 77 A (t' ~ 3 au for Case A, 
and - 7.5 au for Case B) based on ID cmegen models of 77 Car by 
iHillier et ^ (l200lLl200q) . We do not include the inner He^^ zone 
in 77 a’ s wind in the r = 10 a simulations since it is of negligible size 
(r < 1 au), and because 77 B produces essentially zero He^-ionizing 
photons. Due to their smaller domain size, for added accuracy, the 
23 au simulations do include the He^^ structure in 77 a’ s innermost 
wind (r ~ 0.7 au for Case A , and ~ 0.8 au for Case B). 

As discussed in IC14bl our results at apastron appear to rule 
out low values of such as that in Case C. This is supported 
by the work of lM13l and recent observations of 77 Car (T. Gull and 
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M. Corcoran, private communication). Therefore, we focus in this 
paper on simulation Cases A and B, which represent the most likely 
values of t/a’s current mass loss rate. Moreover, little information 
would be obtained by including Case C since for such a low mass 
loss rate t/a already singly-ionizes He throughout the entire com¬ 
putational domain in the r = 1.5 a simulations, and nearly the entire 
domain (out to r ^ 120 au) in the r = 10 a simulations. 

Finally, we note that we use a standard Cartesian coordinate 
system with the origin located at the system centre of mass and 
the orbital plane set in the xy plane, with the major axis along the 
x-axis. The stars orbit counterclockwise if viewed from along the 
-hz-axis. In this system of reference, at periastron, t/a is to the right 
and t/b is to the left. 

2.1 Initial ionization state of the gas 

As in lC14bl. we post-process the SPH simulation output. This ap¬ 
proach works well for the 77 Car system at apastron, when the ion¬ 
ization state reaches an equilibrium value on a time-scale much 
shorter than the orbital time-scale. More importantly, at apastron, 
the t/b wind cavity is an almost axisymmetric cone, allowing 
t/b’s ionizing photons to completely ionize the low-density sec¬ 
ondary wind. In other words, it is reasonable at apastron to assume 
that knowledge of the ionization structure of the preceding SPH 
timestep is not required for accurate RT results. 

In this paper, however, we deal with phases at periastron 
where, due to rapid orbital motion, the interacting wind structures 
evolve on a time-scale comparable to the time necessary for the 
RT simulations to reach a stable ionization state. Moreover, the t/b 
wind cavity at periastron is strongly asymmetrical and distorted. 
Regions that were ionized by t/b in an earlier snapshot may not 
be reachable by t/b’s photons in the next snapshot. To tackle these 
issues, we use a slightly different approach than that in lC14lJ. 

Using collisional ionization equilibrium, the results of lC14bL 
and estimates for the recombination time of He^ as a function of 
gas density and temperature, we are able to better constrain the 
initial ionization state of the gas for use as an initial condition 
in our periastron simplex simulations. We know that He^^ in the 
cooler {T ^10^ K) unshocked gas in the t/b wind on the apastron 
side of the system is completely ph otoion ized by t/b during most 
of the binary orbit around apastron (IC14bh . Observations with the 
Hubble Space TelescopelSp^iCQ Telescope Imaging Spectrograph 
(HST/STIS) show that numerous high-ionization lines that require 
t/b’s ionizing flux are present for most of the orbit and do not start 
to fade u ntil phase ^ ^ 0.984 (a pproximately one month before pe¬ 
riastron, iGulletal . 11200 ^ lM13h . Therefore, we can safely assume 
that t/b is able to completely ionize He^^ throughout its unshocked 
wind at least until 0 ^ 0.985 (and at least for the length scales of 
interest here, r < 155 au from the stars, since the HST observa¬ 
tions show ionized structures extending several thousand au from 
the stars). After 0 ^ 0.985, t/b becomes deeply embedded in the 
dense wind of t/a and the WWIR becomes so distorted by orbital 
motion that we can no longer assume that t/b can efficiently ionize 
He^^ throughout its entire cool unshocked wind. If the He in the 
low-density (n < 10^ cm“^) unshocked t/b gas is no longer pho¬ 
toionized by the star, then, at T = 10^ K, He^ should recombine on 
a time-scale tree ~ 42.4 days (or longer). Thus, around periastron, 
for gas with T ^ 10^ K and 10^ cm“^, the recombination time- 
scale for He^ is about an order of magnitude larger than the orbital 
time-scale. 

Any gas that recombines very slowly compared to the orbital 
time-scale around periastron poses a potential problem for the accu¬ 


racy of our RT simulations, mainly because we are post-proeessing 
the SPH simulation output. The SPH output contains no informa¬ 
tion about the previous ionization state of the gas before we start 
our SIMPLEX simulations. However, any gas with a low enough re¬ 
combination rate should remain ionized from one snapshot to the 
next, even though t/b photons no longer reach that material- There¬ 
fore, we should preserve the ionization state of low-recombination- 
rate gas in the initial condition of our RT simulations. 

The time difference between 0 = 0.985 and periastron is 
^ 30 days. The recombination time for He^ gas with T ^ 10"^ K and 
a ^ 10^ cm“^ is about two weeks longer than this. Lower-density 
gas at r ^ 10^ K will have even longer recombination times. If we 
assume that the t/b ionizing flux becomes trapped by t/a’s dense 
wind and the WWIR at 0 0.985, then at periastron, any such 

low-recombination-rate gas that was initially photoionized by t/b 
at 0 ^ 0.985 should still be ionized. 

Therefore, to mimic the fact that the low-density, cool t/b 
wind on the apastron side of the system should remain ionized at 
periastron, we set the He fractions in the t/b wind material with 
n < 10^ cm“^ and T < 4x 10"^ K to be ~ 97% He^. This should 
generate an initial condition closer to the previous ionization state 
of the photoionized t/b gas, since for gas with a < 10^ cm“^ and 
r < 4x lO'^ K, the He^ recombination time-scale is ^ 4 days or 
longer. The use of collisional ionization equilibrium as an initial 
condition for the remaining gas remains a valid approximation even 
around periastron, since for the pre- and post-shock t/a wind, the 
high densities allow the gas to recombine almost instantaneously. 
Hot (T > 4x 10^ K) gas in the post-shock t/b wind is collisionally 
ionized to He^ or higher. 

Due to the shorter dynamical time-scales around periastron, 
which alter the gas distribution from snapshot to snapshot, we must 
also run simplex for a m uch sh orter simulation time, compared to 
that used at apastron in IC14bl . In our simulation testing we find 
that, setting the initial condition of the gas as described above, in 
Case A, the ionization state reaches an equilibrium value in less 
than one day. In Case B, the gas reaches ionization equilibrium 
in ~ 2-3 days. Since significant changes to the gas density and 
temperature distributions around periastron also occur on a time- 
scale of ~ 2-3 days, we use a total RT simulation time of 2.5 days 
and a simulation time-step of ~ 5 minutes in order to achieve an 
accurate RT calculation of the ionization volumes and fractions. At 
these time-scales, the gas distribution, even around periastron, can 
be assumed to be roughly constant. 


3 RESULTS 

To help the reader better understand the RT results, we first briefly 
describe the density and temperature structure in the xy, xz, and yz 
pla nes. F or more details on the hydrodynamical simulation results, 
see iMll Figs. [2 and El show the number density (top row) and 
temperature (bottom row) for the 155 au and 23 au Case A simu¬ 
lations, respectively (the green circles in Fig. [T] mark the domain 
size of the 23 au simulations). As described in 1^ approach¬ 
ing periastron, the orbital speed of t/a relative to t/b inereases to 
^ 360 kms“^, close to t/a’s wind speed (420 kms“^), highly distort¬ 
ing the WWIR. The small stellar separation at periastron (~ 1.5 au) 
also prevents t/b’s wind from reaching its terminal speed before 


^ Of course, this depends on how long it has been since the supply of ion¬ 
izing photons has been cut off. 
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Figure 1. Slices in the orbital xy (left column), xz (middle column) and yz (right column) planes through the larger (10a) 3D simulation volume for the Case A 
simulation of 77 Car at periastron. Rows show the number density (top, log scale, cgs units) and temperature (bottom, log scale, K). The length scale is shown 
under the bottom left panel. In the first column (i.e. the orbital plane) t/a is to the right and t/b is to the left. The green circle marks the boundary of the smaller 
(r = 1.5a = 23 au) domain simulations (see Fig.O. 


colliding with ha’s w ind. Combined with the effects of radiative in¬ 
hibition by t/a IMIJ), this leads to a significant slowing of t/b ’s pre¬ 
shock wind, which alters the wind momentum ratio at the WWIR. 
This decreases the WWIR opening angle and moves the WWIR 
apex closer to t/b . This behaviour occurs in the Case B simulations 
in Figs. [ 3 ] and m as well, although lowering further increases 
the WWIR opening angle and moves the WWIR apex closer to t/a- 

The above effects create a high asymmetry between the trail¬ 
ing and leading arms of the WWIR at periastron, as t/b becomes 
deeply embedded in t/a’s dense inner wind. The only direction t/b 
can effectively drive its wind is away from t/a (top left panel of 
Figs. [ 2 ] and ID). This has a strong impact on the location of the 
shock-heated gas. As shown in the bottom left panel of Figs.[2]and 
m the post-shoek secondary wind is heated to temperatures higher 
than 10^ K only where the t/b wind collides with the highly dis¬ 
torted leading arm of the WWIR (in red). The trailing wind of t/b 
is unable to collide with t/a’s downstream wind, and so there is no 


hot shocked gas there. This effect is also visible in the bottom left 
panel of Figs. [T] and [3] However, in the larger domain simulations, 
some residual hot gas remains at the outermost right edge of the 
simulation domain in the xy and xz plane slices. This is the expand¬ 
ing, adiabatically-cooling remnant of the trailing arm of the WWIR 
from just before periastron. 

Looking at the xz plane (middle eolumn of Figs. [DS, dur¬ 
ing periastron, t/b passes behind t/a, allowing, for a short period, 
the primary wind to expand in the apastron (-hx) direction. During 
this time, the apex of the low-density cavity on the apastron side of 
the system fills with dense primary wind. The thickness and den¬ 
sity of this inner primary wind region increases with since 
higher move the apex of the WWIR farther from t/a, allowing 
the primary wind to fill a larger volume in the apastron direction. 
Slices in the xz plane also sample mostly cold (10"^ K) gas from 
both winds, although the t/b wind cavity in the Case B simulations 
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Figure 2. Same as Fig.^ but for the smaller (i.e. 1.5a) simulation domain. The black dots in the centre of the top left panel mark the positions of the stars. 


contains warmer materi al due to the larger WWIR opening angle 
and less oblique shocks (lM13h . 

Slices in the _yz plane (right-hand column of Figs.[T}|4]) simply 
show that at periastron, the leading arm of the WWIR has passed 
through the yz plane. There are now clear differences in both the 
density and temperature along the y axis. The -y side of the sys¬ 
tem still consists of unshocked t/a wind, but to the -hy side there 
is a small eavity filled with hot, shocked t/b wind that is bordered 
by dense post-shock t/a wind. This cavity is larger in the Case B 
simulations due to that simulation’s larger WWIR opening angle. 


3.1 He ionization structure and influence of 

Figs. [5] and [6| present the He ionization structure in the orbital xy 
plane for Cases A and B, respectively. In all figures, columns il¬ 
lustrate, left to right, the fractions of He^^, He^, and He^^. Rows 
display the 23 au domain without the inner ionization volume due 
to t/a (top), the same domain with the inner He^ and He^^ ion¬ 
ization volumes (middle, filled yellow and blue circles in the left¬ 


most panel, respectively), and the larger (155 au) simulation do¬ 
main with the inner He^ ionization volume due to t/a included (bot¬ 
tom). The green contour highlights the boundary between the un¬ 
perturbed pre-shock t/a wind and the density-enhanced post-shock 
t/a wind. Due to the complex situation at periastron, the green con¬ 
tour is meant more as a guide to better understand the results than 
a preeise position fo r the pr e-shock t/a wind. 

As discussed in lC14bl. t/b should not have many He^-ionizing 
photons and, therefore, He^^ is produced principally through colli- 
sional ionization. As expected, neglecting the inner He^^ region in 
t/a’s wind, He^^ is only present in the hottest {T > 10^ K) regions 
of the system, i.e. the shock-heated gas in the leading arm of the 
WWIR, and the adiabatically-cooling hot gas in the remnant of the 
trailing arm (from before periastron p assage ). 

Similar to the case at apastron (IC14bh . t/b ionizing photons 
create different He^ structures depending on the value of . In 
Case A, the He^^-ionizing photons are only able to ionize a thin 
layer within the leading arm of the WWIR, and a rather small vol¬ 
ume of the diffuse t/a wind that is expanding into the t/b wind cavity 
located near the remnant of the WWIR’s trailing arm (first row of 
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Figure 3. Same as Fig.^ but for the Case B simulation. 


Fig.EJ. The dense primary wind expanding into the t/b wind cavity 
is able to stop t/b’s He^^-ionizing photons. The He^ present in the 
cooler, unshocked t/b wind on the apastron side of the system is the 
previously photoionized gas that is slowly recombining. The inclu¬ 
sion of the inner He^ and He^^ regions due to t/a have basically no 
effect on the larger-scale He ionization structures in Case A (com¬ 
pare the top and middle rows of Fig. [5]). The only noticeable differ¬ 
ence is that the He^^-ionizing photons are able to penetrate through 
the leading arm of the WWIR and further ionize a very small area 
inside the inner 3 au t/a He^ ionization volume (the red within the 
yellow circle, close to the centre in the -y direction in the first panel 
of the middle row). For Case A, the dense outer primary wind (i.e. 
the area beyond the green contour) is almost entirely He^^. 

The larger WWIR opening angle in Case B, due to the lower 
Mr]p,^ allows the He^^-ionizing photons to escape the inner core and 
singly-ionize He in a portion of the the expanding primary and sec¬ 
ondary winds on the apastron side of the system (right-hand side of 
the panels in Fig. n. Contrary to Case A, introducing the regions 
of ionized He in t/a’s inner wind changes the extent of the ion¬ 
ization front in both the primary and secondary winds. The larger 


(7.5 au) He^ zone in Case B allows ionizing photons to penetrate 
deeper into the receding primary wind, producing a region of He^ 
(although the fraction of He^ remains quite low, see the red region 
in the middle column in the middle and bottom rows of Fig. [6]). The 
He^ on the apastron side of the system is a combination of ongoing 
photoionization and previously ionized gas which is slowly recom¬ 
bining. However, as in Case A, the dense primary wind is mostly 
composed of He^^ and very effectively traps t/b’s He^^-ionizing 
photons (bottom row of Fig. H. 


4 DISCUSSION 

From the discussions in Sections|2]and[3l it is evident that the inter¬ 
pretation of the results at periastron is much more difficult then at 
apastron, especially in terms of trying to understand h ow the results 
relate to observations of the various He line profiles (iNielsen et al.l 
I 2 OO 7 I : IPamineli et ^l2008ah . Nevertheless, we can use our results 
to try to better constrain the regions where the He emission and 
absorption features are generated. We focus on the Case A Simula- 
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Figure 4. Same as Fig.|2] but for the Case B simulation. 


tion results, since as discussed in lM13l : [ci4bl. and above, they most 
likely better represent the current value of Mr ]a • We note that spec¬ 
troscopic phase, defined b y e.g. the disappearanc e of the narrow 
component of the Hei lines (iDamineli et al.ll20083 . or alternatively, 
the m inimum in the observed X-ray light curve, ICorcoran et al.l 
ImI, may differ from true orbital phase by up to one month. 

Our results show that, at periastron, t/b is able to singly-ionize 
He in only a very thin layer of the post-shock t/a wind along the 
WWIR leading arm, and a very small region close to the stars. The 
denser t/a wind, expanding into the secondary wind cavity, is ef¬ 
fectively able to stop the He^^-ionizing photons. Therefore, the 
lower-density receding secondary wind on the apastron side and the 
adiabatically-cooling trailing arm of the WWIR from before peri¬ 
astron are shielded from ionizing photons from t/b . We do not ex¬ 
pect this material to be neutral due to the long recombination time- 
scales for gas at these densities {n < 10^ cm“^) and temperatures 
(> 10^ K). The disappearance of the shock-heated region in the 
trailing arm of the WWIR, caused by the curvature of the WWIR, 
causes a strong asymmetry in the high temperature gas which is 
now mainly present in the WWIR leading arm. The He in this very 


hot region (T > 10^ K) is collisionally doubly-ionized. With the ex¬ 
ception of the inner He^ and He^^ ionization regions due to t/a’s 
ionizing fiux. He in the extended primary wind is mostly neutral. 
This includes post-shock gas in the remnant of the trailing arm, and 
a significant portion of the extended leading arm, of the WWIR. 

Observations show that the broad Hei lines are blueshifted for 
most of the orbital period, both in emission and absorption. The 
P Cygni absorption starts increasing approximately three weeks 
before spectroscopic phase zero, reaching maximum absorption 
around phase zero, before de creasing to a complete disappear¬ 
ance y 20 days after phase zero (iNielsen et al.ll2007l : lDamineli et alJ 
l2008bh . The radial velocity of the absorption reaches its maximum 
blueshifted value (-610 kms“^) shortly before periastron, and then 
suddenly shifts during periastron down to ~ -250 kms“^. This ab¬ 
sorption must be produced by material between the observer and 
t/a, either in the WWIR or in the pre-shock t/a wind. 

Because the radial velocity variations seen in the Hei 27067 
absorption are reminiscent of those o bserved from a star in a highly 
eccentric orbit. INielsen et akl (l2007h assumed that the Hei absorp¬ 
tion is associated with the pre-shock wind of t/a, ionized by the 
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Figure 5. Slices in the orbital plane through the 3D simplex simulation volume for the Case A simulation. Columns, from left to right, show the computed 
fractions of He^"^, He"^, and He^^ (log scale). Rows show, from top to bottom, the smaller (23 au) simulation domain without the inner ionization volumes due 
to the primary, the smaller (23 au) simulation domain with the He^ and He^^ inner ionization volumes due to the primary included, and the bigger (155 au) 
simulation domain including the He"^ ionization volume due to the primary. In this and Fig.|6] the green line marks the edge of the pre-shock primary wind. 


far-UV flux from T/p. However, this assumption leads to an un¬ 
usual, extreme system mass ratio of ^ ~ 0.1, with a mass for t/b 
of 210 Mo and a mass for t/a of ~ 20 Mq, in strong disagreement 
with the most commonly accepted values. In contrast, based on the 
maximum velocity of the absorption seen in the He i line and the re¬ 
quireme nts of high excitation flux and gas density, IPamineli et al.l 
(l2008bl) conclude that the only plausible locations for the forma¬ 


tion of the broad Hei absorptions are the walls of the WWIR. Our 
RT results favour this second scenario. At periastron, t/b is deep 
within the primary wind and its He^^-ionizing photons are not able 
to reach the pre-shock t/a wind or penetrate beyond the dense layer 
of post-shock 77 a wind. As shown in Fig. |5] (central column, middle 
and bottom rows). He in the pre-shock t/a wind is almost entirely 
neutral. The region of He^ is conflned to the centre of the system. 
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Figure 6. Same as Fig. [5] but for Case B. 


Modeling of rj Car’s X-ray light curve and observed broad, 
extended emission from numerous forbidden emission lines have 
helped constrain the orbital inclination (/), argument of periap- 
sis (tu), and sky position angle of the rj Car binary sy stem 

dOkazaki et alJf2()08] : IParkin et al.ll2009l : lMadura et al.ll20123l . It is 

generally agreed that the rj Car binary is inclined away from the 
observer by roughly ~ 45°, with oj ^ 240°-285° and PA^ 302°- 
327°, which places t/b at apastron on the observer’s side of the 
system and implies that t/b orbits clockwise on the sky. With this 
orientation of the binary on the sky, at apastron, our line-of-sight 


(LOS) is more perpend icular to the apex of the WWIR surface. 
Based on the results in IC14bl . we propose that the much denser 
WWIR (specihcally, the compressed post-shock t/a wind) should 
dominate the P Cygni absorption at phases around apastron. How¬ 
ever, because our LOS is more perpendicular to the WWIR around 
apastron, the column density of He^ material between us and the 
continuum source t/a is relatively modest, resulting in a moder¬ 
ate amount of absorption. Moreover, because we are viewing the 
WWIR inclined away from us at roughly 45 degrees, material flow¬ 
ing toward us within the wall of post-shock t/a wind will have a 
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maximum LOS velocity of ~ 300 kms”', since the maximum true 
velocity of the gas within the compressed wall of post-shock t/a 
wind will be roughly equal to the terminal velocity of t/a’s wind. 

When approaching periastron, the trailing arm of the WWIR 
sweeps across our LOS. When this occurs, our LOS is parallel to 
and intersecting the WWIR. This results in a significantly increased 
column of He^ between us and t/a, producing an increase in the 
amount of Hei absorption. Because we are now looking directly 
down the WWIR, we are also seeing the fastest material in LOS, 
resulting not only in an increase in the amount of Hei absorption, 
but also an increase in the blue-shifted velocity of this absorption. 
This scenario is in qualitative agreement with both the increase in 
the absorption and its shift toward more blue-shifted velocities of 
the observed Hei P Cygni profiles going into periastron. Just af¬ 
ter periastron, the strong Hei X1Q61 absorption component rapidly 
vanishes. This is also consistent with our results, since we see that 
at periastron, t/b is no longer able to ionize the WWIR or the pri¬ 
mary wind. With the observer on the apastron side of the system, 
the gas in LOS, further away from the central region, cannot be 
ionized by t/b, resulting in effectively no He^ between us and the 
central t/a He^ reg ion, an d thus a lack of Hei absorption. 

The results in lC14hl imply that the Hei emission seen around 
apastron should arise from a combination of the inner t/a He^ zone, 
the pre-shock t/a wind close to the WWIR apex ionized by t/b, and 
the layer of post-shock t/a wind in the WWIR that is also ionized 
by t/b. The different regions, depending on the orbital phase, where 
the Hei emission can form might help expl ain the multiple emis sion 
components observed in the lin e profiles (iNielsen et al.ll200^ . At 
periastron, as discussed in lMlSl. the WWIR apex is inside t/a’s He^ 
region. The fact that t/b is able to ionize an extremely small region 
close to the centre of the system during periastron means that the 
dominant source of Hei emission at periastron is the inner t/a He^ 
region. Based on our models, little to no Hei emission is expected 
from the WWIR at and just after periastron. 

Just before periastron, the Hen A46S6 line intensity increases 
suddenly and then drops sharply to zero, after which it recovers to 
a second peak before declining back to zero ( Steiner & Daminelil 


I 2 OO 4 I : iMartin et “^120061 : iMehner et al.ll201 ll : iTeodoro et al.l 2012h . 

Based on the density and energy required for the formation of this 
line, the most plausible region in which it c an form during pe¬ 
riastron passage is close to the WWIR apex (iMartin et al.l l2006l : 
iTeodoro et al .1 [20 1 3) . Our results support the scenario proposed 
by ImISI in which the Hen 44686 emission, observed at phases 
0 ~ 0.98-1.3, is generated in the inner He^ region of t/a’s dense 
wind, with the necessary He^-ionizing photons coming from the 
radiatively-cooling shocks in the WWIR. Our RT results show that, 
near periastron, t/b ’s ionizing flux is unable to produce a significant, 
extended region of He^ in either the dense pre-shock t/a wind or the 
densest parts of the WWIR. Even though the unshocked t/b wind 
remains in the He ^ state during periastr on, as discussed by e.g. 
IMartin et*^ (l2006h . lTbodoro et ^(l2012h . the much lower-density 
t/b wind is incapable of adequately explaining the observed Hen 
44686 emission around periastron. Thus, the only remaining region 
that can give rise to the observed Hen 44686 emission is the He^ 
zone in t/a’s dense inner wind. 

The sharp drop in Hen 44686 emission seen near periastron 
might be caused by a combination of two effects. First, we note 
that at periastron, the WWIR opening angle decreases consider¬ 
ably due to the much lower wind momentum ratio, caused by a 
large decrease in the pre-shock velocity of t/b’s wind due to the 
decreased orb ital se paration and radiative inhibition effects (see 
Section [3] and lM13h . At periastron, t/b and the narrowed WWIR 


are also behind t/a and its optically-thick wind (relative to the ob¬ 
server). Thus, t]b and the WWIR are both ‘eclipsed’ by t/a and 
its dense wind at periastron. Second, if the WWIR ‘collapses’ at 
peri astron due to the reduced y/B wind speed (se e discussions in 
e.g. lCorcoran et al ]|201Ql . l^rkin et al.ll201ll : lMO . the radiatively- 
cooling WWIR shocks might disappear. If this occurs, the source 
of He^-ionizing photons (i.e. the soft X-rays generated in the ra¬ 
diative shocks) necessary to produce the Hen emission might also 
disappear. There would then be no significant Hen emission at pe¬ 
riastron. We believe that a combination of a wind-eclipse by t/a and 
a ‘collapse’ of the WWIR is currently the best scenario for explain¬ 
ing the observed drop in Hen 44686 emission near periastron. 

The Case B simulation results are similar to those of 
Case A, but with two ma jor differences. The ID cmfgen models 
(iHillier et al.l l200lL l2006h predict a radius for the inner He^ re¬ 
gion, due to t]a ionization, of ~ 7 au, more than a factor of two 
larger than in Case A. Even without this larger ionization region, 
t/b photons are able to ionize the denser primary wind that is ex¬ 
panding into the t/b wind cavity and reach the receding unshocked 
t/b wind on the apastron side of the system. Qualitatively, we ex¬ 
pect these differences to produce variations, compared to Case A, 
both in the strength and velocity of the He emission and absorption 
lines. Moreover, the fact that t/b ionizing photons are still able at 
apastron to reach the periastron side of the system, have important 
implications on the phases at which these observed features occur. 
The increase in WWIR opening angle as drops should also al¬ 
ter the timing of when specific He emission and absorption features 
are observed. 

Finally, we wish to comment on one detail not implemented 
in our 3D simulations that may alter slightly the He ionization 
structure around periastron. Our current models do not include an 
ionizing s ource at the loc a tion of t/a- The reasons for this, dis¬ 
cussed in IClementel et al.l (l2014ah . are related to the complexi¬ 
ties of proper spectroscopic modeling of the dense, extended wind 
photosphere of t/a, which requires inclusion of the many ioniza¬ 
tion states of metals up to and including iron. However, we note 
that during periastron passage, t/b and the WWIR penetrate deep 
within the extended wind of t/a, down into its hotter regions where 
He is ionized to He^. This ‘bore-hole’ effe ct has many reper¬ 
cussions for multiwavelength observaMes (seelMadura & Owockil 
l201Ql.iMadura et al.ll2012a[lGroh and it may affect 

the He ionization structure of the innermost regions of the rj Car 
system as well. This is because the hole created by the WWIR cav¬ 
ity within the innermost t/a wind may provide an escape path for 
higher-energy photons from t/a- These higher energy t/a photons 
may be able to escape through the WWIR cavity and help ionize 
He in portions of the dense arms of the WWIR, and possibly a 
very thin layer of pre-shock t/a wind located beyond the WWIR. 
Such an effect would be limited to around periastron and extremely 
short-lived since the WWIR geometry changes so rapidly. Escaping 
higher-energy photons from the interior of t/a would only be able 
to escape in very specific directions and travel a limited distance 
before the WWIR geometry shifted and shut-off (or changed the 
direction of) the escaping radiation. Thus, the effect would be lim¬ 
ited to the innermost regions. It is also questionable whether such 
changes to the He ionization structure would even be observable, 
since the WWIR and bore-hole would be located on the opposite 
side of t/a and directed away from the observer. More detailed 3D 
radiation-hydrodynamics simulations are necessary to accurately 
determine how much this bore-hole effect may alter the inner He 
ionization structure and affect observations of different He lines. 
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5 SUMMARY AND CONCLUSION 

We presented the ionization structure of helium at periastron in 
the inner regions of the 77 Car binary system due to the hot 
secondary star’s ionizing flux and collisional ionization. In our 
analysis, we considered two different Mr]^^ Case A {Mr ]a - 
8.5x10-4 Mo yr"') and Case 4.8x10-4 Moyr'). We 

created He ionization maps by post-processing 3D SPH simulations 
of the rj Car system with simplex. Below we summarize our most 
important results and conclusions. 


1. At periastron, in the Case A simulation, the dense t/a wind 
expanding into the t/b wind cavity is able to stop t/b’s He^^- 
ionizing photons. Helium is only singly-ionized by t/b in a very 
thin layer along the WWIR leading arm, and in a very small 
region close to the stars. 

2. In Case B, the larger WWIR opening angle and lower allow 

t/b to ionize a region of the lower-density unshocked t/b wind 
expanding in the apastron direction. 

3. The inner He^ volume due to t/a does not produce relevant dif¬ 
ferences in the final ionization structures due to t/b in simulation 
Case A. In Case B, the larger t/a He^ radius allows the t/b ion¬ 
izing photons to reach deeper into the dense t/a wind and create 
a lower ionization fraction in the pre-shock primary wind. 

4. Collisional ionization creates He^^ in the hot shock-heated gas 
in the leading arm of the WWIR and in the adiabatically-cooling 
remnant of the trailing arm formed just before periastron. 

5. Our RT results support a binary orientation in which apastron is 
on the observer’s side of the system, with the trailing arm of the 
WWIR sweeping across our LOS as periastron is approached. 

6. The small He^ region at the centre of the system is the main 
region were Hei emission can form during periastron. There is 
little to no Hei absorption at periastron because t/b’s ionizing 
flux is effectively trapped by t/a’s dense wind, leading to a lack 
of dense He^ gas between us and the continuum source t/a- This 
is consistent with the observed changes in strength and radial 
velocity of the Hei P Cygni lines across periastron. 

7. A lack of any significant He^ in/near the WWIR at periastron, 
due to t/b, implies that the increased Hen 44686 emission that 
is observed around periastron must arise from the inner He^ re¬ 
gion in 77 a’ s dense wind. Our results favour a scenario in which 
a combination of a wind-eclipse by 77 A and a ‘collapse’ of the 
WWIR are responsible for the observed drop in Hen 44686 
emission at periastron. 


Our model He-ionization structures at periastron are in quali¬ 
tative agreement with available observations of 77 Car’s He lines. 
In the future, we plan to generate synthetic He line profiles for 
comparison to available observational data in order to place tighter 
constraints on the orbital, stellar, and wind parameters of 77 Car. 
These synthetic spectra, together with those generated for the spec¬ 
troscopic high state around apastron, will help us better understand 
the numerous observed spectral features generated in the inner 
~ 150 au of the system. Our final aim is to use 3D time-dependent 
radiation-hydrodynamics simulations to more properly follow the 
ionization and recombination of the gas without postprocessing. 
This is particularly important in the rapidly changing conditions 
and highly asymmetric situation at periastron. 
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